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Recently computed terms of orders O(afn^) in the perturbative series for the r decay rate, and 
similar (new) strange quark mass corrections, are used to discuss the validity of various optimization 
schemes. The results are then employed to arrive at improved predictions for the complete terms 
order O(a^) and 0{a s a ) in the massless limit as well as for terms due to the strange quark mass. 
Phenomenological implications are presented. 
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I. INTRODUCTION 

The dependence of the r-decay rate on the strong coupling a s has been used for a determination of a s at lower 
energies, with the most recent results of 0.334 ± 0.007 oxp ± 0.021 th oo and 0.3478 ± 0.009 exp ± 0.019 tho o by the ALEPH 
jl| and OPAL [|| collaborations. After evolution up to higher energies these results agree remarkably well with 
determinations based on the hadronic Z decay rate. In view of the relatively large value of a s (M T ) estimates for the 
yet unknown terms of higher orders play an important role in current determinations of a s at low energies. This is in 
contrast to high energy measurements, where the uncertainty from terms of order of 0.001 to 0.002 |3| is somewhat 
less or at most comparable to the present experimental errors. 

The situation is even more problematic for the determination of the strange quark mass from the Cabbibo suppressed 
r decays. Perturbative QCD corrections affecting the m 2 s term are extremely large and contributions from increasing 
powers of a s are barely decreasing, which casts doubts on our ability to extract a reliable result for m s from this (in 
principle) clean and straightforward measurement ||, ^, ^, |[ |l0[ |ll) . 

Partial results of order a* for the absorptive part of the massless vector and scalar correlators have been obtained 
recently namely terms proportional to rij, where rif denotes the number of massless fermion species. These 
allow to test two popular optimization schemes — based on the principles of "minimal sensitivity" (PMS) and of 
fastest apparent convergence (FAC) |l3|, [l4|, [l5| — which have been used to predict yet uncalulated higher order 
terms (T^ [l7j. 

It will be demonstrated that the predictions of both schemes (coinciding at af) for the coefficient of order nja^ 
are in reasonable agreement with our calculations, which are then used to predict the complete fixed order (FO) and 
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the "contour improved" (CI) (T^, |l9| 0{a 4 s ) contributions to the r decay rate. Employing the four-loop QCD beta 
function in combination with improved terms a rough estimate even for 0{oi s ) terms can be obtained. The results 
lead to fairly stable values for a s consistent with current analysises. 

Unfortunately, no essential decrease of the difference between the central values a s (M T ) as obtained with FI and 
CI approaches is observed after including order aj. and corrections, assuming for the moment that these estimates 
are indeed correct. On the other hand, the theoretical uncertainty assigned to a s {M T ) within each method according 
to standard techniques does decrease significantly. 

The implication of this approach for the extraction of m s from Cabbibo suppressed decays will be investigated 
along the same lines. New results will be presented for the terms of order UfOt^ m 2 in the total rate. In this case the 
agreement between FAC /PMS predictions and our results is quite encouraging and naturally suggests the use of the 
former as a reliable prediction for the complete a^m 2 term. Following an approach discussed in Jl7| , even a rough 
estimate Tn 2 a^ terms can be obtained from these considerations. 

However, the rapid increase of the coefficients indicates that the inherent uncertainty of the present m s determi- 
nations will not necessarily decrease with inclusion of the higher orders. As we will see, the situation is somewhat 
better for the spin one contribution if considered separately. 



II. GENERALITIES 



We start with the well-known representation 
integral along a circle C of radius \s\ = M 2 



Hi f22, [23, E4| of the tau-lepton hadronic rate as the contour 



R T = 6in 



ds ( s V 



(i) 



Here II W and nJ 2 ! are proper flavour combinations of the polarization operators appearing in the decomposition of 
the correlators of vector and axial vector currents of light quarks 



n^( gi m iimj ,m, ftas ) =i dxe^(T[j^(x)U^f)H0)}) 



9^l v/A (q 2 ) + q , q ^ y/A { q 2 ) 



(2) 



with m 2 = J2f=uds m2 f an d J^lif = qiJ^(l5)<lj- The two (generically different) quarks with masses rrii and mj are 
denoted by % and qj respectively. 

For the case of the r-lepton the relevant combinations of quark flavours are ij = ud and ij — us. The polarization 
functions Hy/A' I = 1, 2 are conveniently represented in the form (Q 2 = —q 2 ) 



(Q 2 ) (l - 2) rtl v/A (q 2 ) = 



JLttM (V- a ) + — V Q- D U ll] m 2 a) 

^ D>2 ^ 



(3) 



Here the first term on the rhs corresponds to the massless limit while the first term in the sum stands for quadratic 
mass corrections. We neglect the masses of u and d quarks. Therefore in perturbative QCD Ily = II^ and we will 
often omit the subscript V/A in the following. Current conservation implies: = IIq . 



The full tau-lepton hadron rate R T can be presented as a sum of spin 1 and spin parts, viz. 

2 

A/? 



R w =6i7T I A ( i _ _!_ 



f + 2^jn«( s ) + nW(o)A< 



v?'." 1 = ()/ti / — 

S | =M 2 M 2 



n(°)( s ) -nM(o)/s 



(4) 



where 



nd) = -nW/q 2 , n(°) = + n^/q 2 



(•5) 



and the contribution of the singularity at the origin (proportional to IlM(O)) has to be included. A nonvanishing 
value of nW (0) is a nonperturbative constant. 
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On the other hand, the unknown constant drops out if one considers moments 

s ° da ( s \ k ( s \ l dR { T 10) 



"w-w. < 6 » 



with k > 0, I > 1. (Note that the moments introduced in |19 are related to ours as R^ 1 = Rf' k ' 1 + R^ k ' 1 .) 



The decay rate R T may be expressed as the sum of different contributions corresponding to Cabibbo suppressed or 
allowed decay modes, vector or axial vector contributions and the mass dimension of the corrections 



R T — R T ,v + Rt.a + Rt.s (7) 



with 



Rv — -\V u d\ 2 


1 - 


f $o- 


f E $V,ud,D I 








D=2,4,... / 


Ra= l\V ud \ 2 


1 - 


f S Q - 


f E £4,ud,23 








£>=2,4,... / 


Rs = 3|K S | 2 J 


H 


-<5 4 


E &us,D J ■ 








d=2,4,... y 


dimension of the fractional corrections Sy/A,ij,D, 



(8) 



vector and the axial vector contributions: 8%j d = (^v,y,D + ^A.ij,D)/2. If a decomposition into different spin/parity 
contributions is made or a particular pattern of moments is considered then we will use the corresponding obvious 
generalization of (||). For instance, 

4T = ^ik s i 2 [i+5o w + E 0) 

y d=2,4,... y 

and 



4T = ik, s i 2 ( E do) 

Thus, in our notation we have the relation 

The integral in eq. (Q) is, obviously, insensitive, to the (^-independent terms in the polarization functions rig 1 ' and 
n£ ] . This means that without loss of generality we may deal with the corresponding (Adler) /^-functions, viz. 

An important property of the functions , and n^ 2 ' is their scale independence, which implies that they are 
directly related to measurements. 

The Adler functions and have been calculated with C(a^) accuracy, the polarization function IT^ 2 ', however, 
to 0{ot 2 s ) only (see and references therein). 

The (apparent) convergence of the perturbative series for is acceptable, the one for and II^ 2 ' is at best 
marginal. This has led to significant theoretical uncertainties in extracting a s and to a fairly unstable behaviour in 
extraction of m s from r-decays. 

To improve the situation, we have computed the two leading terms in the large rif expansion of the next order, 

i.e. terms of order O(n^af), 0{n 2 a A s ) to D 1] , D [ z ] and 0(n 2 f a 3 s ), 0(n/aig) to n^ 2 '. Our results are described in the 
following section. 




4 



III. FIXED ORDER RESULTS IN a 3 a AND a* 

Using the technique described in jl^, |2^, |2(| and the parallel version of FORM |27|, |2^| , the leading and subleading 
(in rif) terms of the next order in the perturbative series for and have been obtained in the standard 

MS renormalization scheme fl29], |3CT| : 



Dk 



+ < 



11 2 ' 
"12 + 3 C3 



151 19 
162 ~ 27' 

87029 
+ - 



1103 



288 



n f + 

7847 
216 



, 275 

C3 H — 5- Cs 
6 



■24 ""Cs 



262 



25 



6131 203 5 
5832 + 324 ^ + 18 



1045381 
15552 



40655 



Cs 



260 



864 ,J 27 
a s + a 2 s {-0.1153 n f 



Cs 



7 2 +rf [1]4 T7 f +J [1]4 
V T "0,1 n f + "0,0 



1.986} 



a 3 {0.08621 n 2 , - 4.216 n/ + 18.24} 



-0.01009 



1.875 n 2 f + d, 



Vl "/ + "0,0 / > 



»/ 



(12) 



D 



[i] 



8671 



11664 27 



11 2 

y + 3 Cs 



13 r 



5185 39 



Cs 



2 

3257 



1265 



81 

396781 461 



Cs 



n/ 



144 

44273 
972 + ~8T 
2641517 131275 



Cs 



5184 



559872 1296 



1 C + 5 C 
48 C4+ 18 Cs 



216 



n 3 f + 



12845 
36 



61913567 59 



1119744 54 



352549 67 22859 
7776 ^ 3 + 9f> + 3888 Cs 



^2,0 



1.667 a s + a 2 s {-0.5736n/ + 12.57} 
a 3 s {0.1646 - 14.31 n/ + 149.} 



a 4 { 



-0.01563 n 3 f + 6.067 n 2 f 



(13) 



5 



n 



-Ami 



7 
3 



+ at 



( 




+ al[ 




' 2131 


+ 


11664 




3997 


486 



25 _ 2 
"24 9 43 



19 



81 



15331 359 520 . 



68135 52 
27 



1944 



3875 
243 ' 



/2 
S3 



J2]3 
^2,0 



-4m 2 1 + 2.333 a s + a 2 s {-1.309 n f + 23.51} 



|o.4647n^ - 32.08 n/ + A^cf} 



(14) 



Here we have used a s 



■,m s 



m s (Q 2 ) and set the normalization scale /i 2 = Q 2 ; results for generic values 



of can be easily recovered with the standard renormalization group techniques. The result for the a 4 terms in D 



[i] 



has been already presented in |12| , the coefficients of the a 4 and terms in in and Ili^ respectively are new. 



IV. IMPLICATIONS FOR THE a 4 PREDICTIONS AND PHENOMENOLOGICAL ANALYSIS 



A. Massless case 



FAC (Fastest Apparent Convergence) and PMS (Principle of Minimal Sensitivity) methods are both based eventually 
on the concept of scheme-invariant properties and the idea of the choice of an "optimal" scheme to provide better 
convergence of the resulting perturbative series. For both methods the optimal scheme depends on the physical 
observable we are dealing with. With FAC it should be a scheme which minimizes (set to zero by construction) all the 
terms of order a 2 and higher, while the PMS scheme is fixed by the requirement that the perturbative expansion for 
the observable is as insensitive as possible to a change in the scheme. The assumption that a renormalization scheme 
is in a sense optimal sets certain constraints on not yet computed higher order corrections in any other scheme. These 
constraints can be used to "predict" (at least roughly) the magnitude of these corrections. 

For the function the result is known since long from Ref. jl6| (see Table 1, column 5). From the three entries 
corresponding to nf — 3,4 and 5 one easily restores the FAC/PMS prediction for the rif dependence of the a A s term 

in Dg 1 ^ (the term of order rvj was fixed to its computed value): 

4 1 ' 4 = 127.6 - 44.2 ti/ + 3.64 n) ~ 0.0100928 rfd 1 ^ 4 = 127.58 - 44.211 n f + 3.6439 n) - 0.0100928 n). (15) 
(Note that FAC and PMS predictions happen to coincide for the term.) 



n f 


jcxact 


dl AC 


d p 3 MS 


^FAC/PMS 


dl AC 


3 


6.371 


5.604 


6.39 


27 ± 16 


145 ± 100 


4 


2.758 


4.671 


5.26 


8 ±28 


40 ± 160 


5 


-0.68 


3.762 


4.16 


-8 ±44 


-3 ± 230 



TABLE I: Estimates for the coefficients ds — dp 1 ' 3 and d± — d^ 4 in the function Dq based on FAC and PMS optimizations. 
The estimate for of rig 1 ' 4 employs the exact value of dj, 1 ' 3 . The last column contains the FAC predictions for the coefficient 
df AC which was obtained assuming the value for dp 1 ' 4 as given in the fifth column; the corresponding uncertainties have been 
estimated as described in the text. 

[113 

It is interesting to compare the FAC/PMS predictions for the n/ dependence of coefficient d with the exact 
result given in eq. [l^. The results of both estimates are 

d [ o ]3 (FAC) = 8.54- 1.013 n f + 0.0116 n^, (16) 
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d [ o ]3 (PMS) =9.93- 1.23 n f + 0.0125 n). (17) 

The comparison of the complete a s s and partial a A s results with FAC and PMS estimates leads to the following 
observations: 

• Starting from a,, the leading in a s and nf terms of order a^nj and a^nj are numerically quite small (at least 
for n/ < 6) and, thus, should have a negligibly small influence on the coefficients of the a s expansion. On 
the other hand, the term subleading in nj, say, of order a^rif is comparable in size with the term of order 
a^n®. Similarly, the a^nj term is significantly smaller than the o^Ti/ one, whereas the a 4 n/ and a A s n°f terms 
are expected to be of similar magnitude. 

• In general FAC/PMS methods correctly reproduce sign and order of magnitude of the higher order coefficients. 
The agreement is getting better when the coefficients happen to be large, as is the case for the quadratic mass 
corrections (see below). 

• Taken separately, the FAC/PMS estimates of the coefficients of the n/ expansion could deviate rather strongly 
from the true result. However, for a given n/, the deviation in the predicted value of the full 0(a") term tends 
to be significantly smaller than what could be expected from summing individual terms of the n/ expansion. In 
addition, for the particular point n/ = 3 very good agreement is observed. 

To illustrate this feature, let us consider the worst case: the FAC/PMS prediction for the o? s nf term in the 
function dffl. Here the ratio of the exact result relative to the predicted one is quite large (about 4). Without 
knowledge of the or s nPf contribution one would expect that the uncertainty of the prediction for the full a 3 
coefficient should be at least around 

CjP-ISi j[1]3i \ 

(«0,1 \exact - «o,l \FAC) nf. 

For nj = 3, 4, 5 this amounts to 9, 12 and 15, which should be compared to the corresponding differences of the 
full order O(a^) (summed over all contributing power of nf) coefficients, viz. 1, 2, 4. This example demonstrates 
that the deviation of FAC/PMS predictions for subleading in nj terms may well serve as conservative estimate 
of the accuracy in the prediction of the complete terms of orders and a*. 



0,2 71 f 



These observations motivate the assumption that the prediction for the coefficient d (= d + d \ nf + 
d^n 3 f) should also be correct within 

, fj[l]4i ,[1]4, s 2 

±(,"0 \exact — «o \FAC) n f 

Thus, in our phcnomcnological analysis of the r-lepton decays the estimate 

4 1]4 |n / =3-27±16 (18) 

will be used. On the basis of this improved estimate and the four loop (3 function |3l|| one may even speculate about 
the term (whose exact evaluation is completely out of reach in the foreseeable future) . Following the discussion of 
Kataev and Starshenko jL7|, one obtains 

4 ll5 | n/ =3 = 145±100, (19) 

not far from the previous estimates of Ref. [l^| . The variation of rff, by ±16 leads to the variation of 4 by ±100. 
For other values of n/, the corresponding predictions can be obtained in the same way. They are listed in Table 1. 

The FAC/PMS prediction (|l5|) for the nf dependence of the coefficient dg 1 ' 4 does not take into account the available 
knowledge of the corresponding rij part. One can easily include this by fitting the FAC/PMS predictions for only 



two values of nf with a linear function of nf. As a result one obtains |38 (we have boxed the predicted coefficients 
in order to separate them clearly from the input) 



4 1]4 (FAC/PMS, n/ = 3,4) = 1 105.7 - 3L8n7|+ 1.875 n) - 0.01009 n), (20) 
4 1]4 (FAC/PMS,n / = 4,5) = I 107.7 - 32.3 n f I + 1.875 n) - 0.01009 n), (21) 



4 1 ' 4 (FAC/PMS, n/ = 3,5) = I 106.4 - 32.0 n f 1+ 1.875 nj - 0.01009 nj. (22) 
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One could now perform a self-consistency check of Table I by predicting, say, d^ 4 for nf = 5 from eq. (^0|). The 
result - (-7.5) - is compared successfully to the value listed in the table, viz. -8. The corresponding predictions from 
eqs. (|l]) and ( p2| ) are also in very good agreement to Table I. 

An instructive example of how knowledge and inclusion of the subleading n / term can improve FAC /PMS predictions 

is provided by the (exactly known) coefficient dp 1 ' 3 . Indeed, assuming the knowledge of cZq 1 ^ 3 and djj^ and using the 
values of <if AC and d^ MS (at n/ = 3) as given by Table 1, one easily arrives at 



"0 
and 



4 1]3 (FAC) = 1 17.48 | - 4.216 n/ + 0.08621 n 2 f (23) 



4 1]3 (PMS) = 1 18.27 | - 4.216 n f + 0.08621 rvj, (24) 



[113 

which should be compared to the exact value d — 18.24. Repeating the same exercise for nf — 4, 5 we get 20.16 
(FAC) and 20.75 (PMS) for n f = 4 as well as 22.69 (FAC) and 23.08 (PMS) for n f = 5 respectively. 

Bearing in mind that in many cases FAC/PMS predictions made for nt — 3 are in better agreement with the exact 
results (see, e.g. Table I and tables from Ref. p7J)j we suggest eq. (20) as the best FAC/PMS prediction for the 

— . [114 — 

constant term and the term linear in rif in the coefficient d . 

Eqs. (|l8|,|l9]) can be used to predict R T in the massless limit first in fixed order perturbation theory (FOPT) 

i?f OPT = (25) 
3(1+ a s + 5.202 a 2 + 26.37 a 3 s + a A s (105 ± 16) + af(138 ± 230 ± 100)) . 

Here a s = a s {M T )/-K. The first uncertainty in the a 5 s term comes from that of dp 1 ' 4 while the second is our estimation 

of the error in the very coefficient fig 1 ' 5 . (Of course, within this approach they are strongly correlated). 

Similarly, we can use "contour improved" (CI) formulae |l8|, [l{| (assuming as reference value a s (M T ) = 0.334 Q|) 
to get 



R? 1 = 3(l + 1.364 a s + 2.54 a 2 s + 9.71 a 3 s + 1.31 a 4 ^ 4 + 0.95 a^ 5 ) (26) 



or, equivalently, 



R^ 1 = 3 (l + 1.364 a s + 2.54 a 2 s + 9.71 a 3 s + a 4 (35 ± 20) + ^(138 ± 95)) . (27) 

Let us compare our new value for the coefficient d ^ 4 in ( |l8| ) with the ones used in extracting a s (M T ) from t data 
by the OPAL @] and ALEPH collaborations, namely 

4 1]4 |„ /=3 = 25 ± 25 (OPAL), 50 ± 50 (ALEPH). (28) 

While OPAL's central value is basically the same as ours, their error bar is somewhat larger. In the case of ALEPH 
both the central value and the uncertainty assigned to it are significantly larger than our numbers. In this connection, 
we would like to stress that eq. ( |l8| ) utilizes completely new non-trivial information given in ( |l^ ) : the subleading term 
in nf term of order af. 

It is of interest to see in detail which accuracy in the determination of a s (M T ) one could achieve assuming eqs. 
(|l8|,[l9|). Let us introduce the quality 8p as follows: 

R T s=o = —^7 f_ °, - = \V u <i\ 2 SewRti (29) 

I (t — > ivv) 

with 

R T = 3(1 + dp + Sew + Snp)- 

The first term here is the parton result, the second stands for pQCD effects. The non-perturbative corrections 
represented by 6np happens to be rather small S NP = —.003 ± 0.003 (see, e.g. [§2|). Here the flavour mixing matrix 
element \V u d\ 2 = 0.9475 ± 0.0016 The factor Sew — 1-0194 is the electroweak correction which collects the large 
logarithmic terms [Q, while Sew = 0.001 is an additive electroweak correction JmJ. Using for definiteness the result 
of ALEPH @ 

R tS=0 = 3.492 ± 0.016, (30) 
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one arrives at 

5 C p P = 0.207 ± 0.007. (31) 



To get a value for a s {M T ) one should simply fit S C p P against R T /3 — 1 as given by eq. ( p5[ ) or by eq. ( p7| ) to get a 
result corresponding FOPT or CIPT (Fixed Order or Contour Improved PT). 

Unfortunately, there is no unique way to assign a theoretical uncertainty 5a s to the obtained value of a s (M T ). In 
the literature one finds several suggestions. Let us consider them in turn. 

1. 5a s is a half of the shift in a s induced by the last fully computed term in the PT (that is by the one of order 
Qs at present). 

2. Sa s is equal to the change in a s caused by varying the normalization point fi around M T , typically within the 
range of 1.1-2.5 GeV. 

The corresponding results read (terms of order and higher in eq. ( |25| ) and in the D-function have been set 
to zero) 

a FOPT (M r ) = 0.345 ± (0.025|0.037), (32) 
a° IPT (M r ) = 0.364 ± (0.012|0.021). (33) 

Here the first (second) value in brackets correspond to the use of the first (second) suggestion for the error 
estimation. After evolution from M T to Mz this corresponds to 



o 



FOPT 

s 

CIPT 



(M z ) = 0.1209 ± (0.0024|0.0037), (34) 
a^ lt " r {M z ) = 0.1229 ± (0.0011|0.0020). (35) 

3. 5a s is equal to the change in a s caused by the uncertainty in the predicted (that is not yet completely known) 
higher order terms in the perturbative series for R T ; 

4. Sa s is a half of the difference in the a s (M T ) as obtained within FOPT and CIPT. This difference comes from 
different handling of higher order terms. 

In order to quantify the error estimates according to 3. and 4. we have shown in Table [n| the results for a s (M T ) 
obtained with various choices for d^ 4 , d^ 5 and /i. The entries with the choices ±100 for the coefficients illustrate the 
large change in a s which would result from a failure of PMS and FAC once higher order terms are included. For the 
plausible values of dp 1 ' 4 and d^ 5 we observe an significant decrease of the \x dependence after inclusion of additional 
terms in the a s series. 

Our final predictions for a F PT (M T ) and ag IPT (M T ) are given in the first column of Table III, together with 



experimental error [|39| and the combined theory uncertainty. The values of theory uncertainties are listed separately 
in columns 3,4 and 5. The corresponding values at the scale of Mz are 

a FOPT (M z ) = 0.1192 ±0.0007 ± 0.002, (36) 
a° 1FT (M z ) = 0.1219 ±0.001 ±0.0006. (37) 



Thus we observe that the total uncertainty based on a combination of not yet calculated higher order terms, /i- 
dependcncc and scheme dependence is reduced, once a A s terms are available. However, the difference between FOPT 
and CIPT results of roughly 0.02 is a remaining, at the moment irreducible uncertainty |40||. 

It is thus of interest to study this difference as a function of m T . In practice, this could be applied to sum rules for 
spectral functions as determined in e + e~ annihilation. 

Therefore, let us consider a hypothetical case of r lepton with the mass equal 3 GeV. Assuming a s (1.77 GeV) = 0.334 
and running this value to 3 GeV via standard 4-loop evolution equation one gets a s (3GeV) = 0.2558 and predicts 

<5p P = 0.1353 

which corresponds to eq. ( p5| ) with the a 4 and a 5 s terms fixed to their central values (see eqs. (|l8|-|l9|) ) . Let us 
now investigate the results for a s and the theory error that would result from 5 C p P = 0.1353 as a starting point 



The corresponding analogs of eqs. (p2H35|) and Table (gl|) are displayed below as eqs. ( (|38|,|39j) ) and Table (IV) 



9 



d [1]4 


at 


d [1]5 
"o 


a 5 3 


27 


0.331 ±0.02 


145 


0.33 ±0.02 


0.357 ± 0.009 


0.354 ± 0.004 


43 


0.329 ± 0.02 


245 


0.324 ±0.01 


0.353 ±0.01 


0.348 ± 0.005 


11 


0.333 ± 0.02 


45 


0.335±? 


0.361 ± 0.008 


0.360 ± 0.002 


100 


0.324 ± 0.02 


100 


0.314 ±0.01 


0.340 ±0.01 


0.338 ± 0.007 


-100 


0.348±? 


-100 


0.349±? 


0.398 ± 0.03 


0.401±? 



TABLE II: The predicted value of a s (M T ) in dependence of chosen values for the coefficients d^ 4 , S , . The second and the 
forth columns differ in the the number of terms in the perturbative series included. The upper value of a s is the one predicted 
within FOPT, the lower corresponds to CIPT. The uncertainty in the value of a s corresponds to changing the normalization 
point n as follows /x 2 /M 2 = 0.4 — 2. The entry with question mark means that an equation for a s (/i) does not have a solution 
for some value of \i within the interval. 



Method 


a 3 (M T ) 


A^ xp 


Afj, 


AdJ 1 ' 4 


Adj, 1 ' 5 


FOPT 


0.330 ± 0.006 ± 0.02 


0.006 


0.019 


0.0045 


0.0011 


CIPT 


0.354 ± 0.009 ± 0.006 


0.009 


0.0036 


0.0042 


0.0019 



TABLE III: The value of a s (M T ) obtained with 8 e p P , dj, 1 ' 4 and df, 1 ' 5 fixed to their central values according to eqs. ( ^l|Jl^ , |l9| ) 
together with corresponding errors. 

correspondingly. The difference between FOPT and CIPT decreases significantly, and this remains true even after 
extrapolating to a s (Mz)- 

0.263 ±(0.013|0.014), (38) 

0.265 ± (0.005|0.008), (39) 

0.1198 ± (0.002|0.003), (40) 

0.1203 ± (0.009|0.0016). (41) 

The same coefficients d^ 4 and d\^ 5 can be used to predict |l(| [l7j (the non-singlet part of) corrections of orders 
and a 5 s to the i?-ratio in Z decays: 

R(n f = 5) = l + a s + 1.409 a% - 12.77 a 3 s ± (-97 ± 44) a 4 ± (76 ± 230) a 5 s . 
It is also of interest to display R(s) for n/ = 3 and 4, which is accessible at e + e~ colliders at lower energies. 

R(n f = 4) = 1 + a s + 1.525 a 2 s - 11.52 a 3 s ± (-112 ± 30) a 4 ± (-245 ± 160)a^. 

R(n f = 3) = 1 + a s + 1.640 a 2 s - 10.28 a 3 s + (-129 ± 16) a 4 ± (-635 ± 100)a^. 

Our results are close to those of Jig, [llj, they employ, however, the additional information from Jl^, |3l| . 

These formulae demonstrate rather good convergency for nj — 5 and a reasonably good one for rif — 3 and 4 if our 

predictions for the coefficients dg 1 ' 4 , dj^ 5 deviate from the true values within the assumed error margins. 



af OPT (3GeV) = 
«P PT (3GeV = 
a FOPT (M z ) = 
a? IPT (M z ) = 



10 



d [1]4 


a 4 


d [1]5 
"o 


a 5 3 


27 


0.256 ± 0.007 


145 


0.256 ± 0.003 


0.262 ± 0.004 


0.261 ±0.002 


43 


0.256 ± 0.007 


245 


0.254 ± 0.003 


0.26 ± 0.004 


0.259 ± 0.002 


11 


0.257 ±0.006 


45 


0.258 ±0.005 


0.264 ± 0.003 


0.263 ±0.001 


100 


0.253 ± 0.008 


100 


0.249 ± 0.004 


0.255 ± 0.008 


0.254 ± 0.002 


-100 


0.264 ±0.01 


-100 


0.277± ? 


0.279 ± 0.02 


0.28± ? 



TABLE IV: The predicted value of a s {M T ) in dependence of chosen values for the coefficients d^ 1 ' 4 , df, 1 ' 5 for a hypothetical 
case of M T — 3 GeV. The third and the forth columns differ in the the number of terms in the perturbative series included. 
The upper value of a s is the one predicted within FOPT, the lower corresponds CIPT. The uncertainty in the value of a 3 
corresponds to changing the normalization point fi as follows /i 2 /M 2 = 0.4 — 2. The entry with question mark means that an 
equation for a s (M fl ) does not have a solution for some value of /x within the interval. 



B. Quadratic mass corrections 

Let us first discuss the function D 2 > . The FAC/PMS predictions can be easily obtained following p7| [; they are 
listed in Table 0. 



n f 


icxact 
"3 


4 AC 


4 MS 


^FAC/PMS 


3 


107.5 


89.86 


91.17 


1200 ± 400 


4 


94.37 


79.13 


80.11 


950 ± 300 


5 


81.54 


68.66 


69.33 


750 ± 200 



TABLE V: Estimations of the coefficients d% = d\^ 3 and dn = dj 1 ' 4 in the function based on the FAC and PMS optimizations. 
The estimation of dj 1 ' 4 employs the exact value of dj 1 ' 3 . 
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We again restore the n/ dependence of the coefficient d l 2 1 as predicted by FAC/PMS: 

4 1]4 (FAC/PMS) = 1931.44- 281.956 n f + 9.0294 n) - 0.0156289 n) (42) 

as well as that of d^ 3 - 

123.654 - 11.6638 n f + 0.133293 nj, (43) 



125.975 - 12.0028?!/ + 0.134769 nj (44) 



148.978- 14.3097n/±0.16463n^. 

The comparison of estimates and exact results reveals a picture qualitatively similar to the massless case but with 
some modifications. A few important observations are in order. 

[113 

1. All three terms of the nf expansion of d 2 are successfully predicted within about 20 % accuracy. 

[113 

2. Unlike the massless case the agreement between FAC/PMS predictions for the coefficient d 2 for n/ = 3,4,5 
and the corresponding exact numbers is within the range of 15%-20%. On the other hand, the estimation of accuracy 



4 1]3 (FAC) = 
4 1]3 (PMS) = 

to be compared with 

d 2 ^ 3 ( exact) = 
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of the fixed n/ predictions obtained exclusively from the knowledge of the subleading contribution of O(o^n/) is 
of the right order of magnitude but somewhat less. All this is probably a consequence of a significantly larger rif 
independent contribution. 

3. At O(a^) the exact result for the full coefficient 



d [1]4 - d [1]4 

u 2 — "2,0 



4,1 n f 



^2,2 n / 



^2,3 n f 



is unknown, apart from its leading and subleading terms in u/. The prediction for the subleading coefficient d 2 ^2 — 9.0 
is by 50% larger than the exact value 6.07. The predicted values for d^Jf an d ^i 4 are ver V large. In view of this 
largeness, the estimate (d^Uxact — d^lFAc) n / (— 30, 50, 80 for u/ = 3,4,5 respectively) looks somewhat too 
optimistic. Therefore we assign a conservative 30% uncertainty to the fixed nt predictions listed in the fifth column 
of the Table [v|. 



n } 
3 


fc 3 FAC 
199.1 


fc 3 PMS 
200 ± 60 


7 PMS 

2200 ± 1500 


4 


171.2 


170 ± 50 


1800 ±1100 


5 


144.7 


145 ± 40 


1400 ± 900 



TABLE VI: Estimates of the coefficients kz = k 



[2]:) 



The estimate of fcl 1 ' 4 employs the predicted value for kz 



and &4 = fci, 2 ' 4 in the function Ilj 2 ' based on the FAC and PMS optimizations, 
the corresponding uncertainties include only the ones induced by kz. 



At last, we repeat the analysis for the function Lit, . The results of FAC/PMS optimization methods are given 

i™ s from the table for nf = 3,4, 5 we reconstruct the corresponding full nf 



in Table VI. Using the values of k% 



dependence: 



[1]3 



294.472- 33.2429 n f 



0.696598 n). 



The comparison with the known terms of order n/ and n 2 (—32. 0843 n/, 0.464663 nf 2 ) demonstrates a remarkably 
good agreement for the subleading nf contribution. The 50% error in the predicted value of the n 2 contribution looks 
natural as the corresponding coefficient is small. Following the same line of reasoning as above we have assigned a 
30% uncertainty to the 0(a 3 ) fixed nj result. 

To get a general idea about the size of the contribution to eq. ( [li]) we used FAC/PMS and the predicted a 3 , 
coefficient. The results are listed in the forth column of Table VI. 

Let us now consider the effect of a 3 and corrections on the determination of the strange quark massEJ. The mass 



correction to the R T depend on both functions D^ 1 ' and Ili^. Let us use the central ALEPH value of a s (M T ) = 0.334 
when estimating the size of the perturbative corrections. For fixed order one finds the mass correction to the total 
rate: 



[2] 



us,2 



= -8— + 5.33 a s + 46.0 + 284 
M T 2 

+ 0.75 a 3 s k l 2 2]3 + a^(723. + 0.25 d[ 1]i + 9.84 fcf 13 + 0.75 k[ 2]4 )) 

777, 

= -8—|(l. + 0.567 + 0.520 + 0.521 ± 0.05 ± 0.593) 

Mi 



-8-^(3.2 ±0.6), 



(45) 



where in the last line we have assumed the (maximal!) value of the 0(a A s ) term as an estimate of the theoretical 
uncertainty (this convention will be used also below). 
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For the "contour improved" series one obtains 

Ml 



$u°s2 = -8^|(1.44 + 3.65 a s + 30.9 + 72.2 + 1.18 a^' 21,3 



+ a*(0.678 4 1]4 + l-06 4 2 ' 4 )) 



-8—^(1.44 + 0.389 + 0.349 + 0.371 ± 0.09 + 0.403) 

Mi 



= -8^|(2.95 ±0.4). (46) 

Now we consider the contributions of spin 1 and spin separately. The lowest moments (L = 0) of the spin- 
dependent functions depend on a nonperturbativc quantity and, thus, can not be treated perturbatively in principle 

0- 

For the spin one part and for I) = (0, 1) we find 

2 

& = - 5 ^|( L + 4 ' 83a + 35.7^ + 276. at + a 4 s (U50 + 
mi 



-5^(1. + 0.514 + 0.404 + 0.331 + 0.326) 

Mr 



2 

mi 



and 



-5^(2.58 ±0.33) (47) 



= -5^(1.37 + 2.55a + 16.1 + 135 a 3 + 0.895 a^4 1]4 ) 



-5—^(1.37 + 0.271 + 0.182 + 0.163 + 0.137) 

Mi 



2 

mi 



= -5^(2.12 ±0.14). (48) 

Note that spin 1 contribution is determined by the component III 1 ! alone and is known up to third order. Clearly, this 
series is decreasing in a reasonable way (comparable to the behaviour of 2 ) and, at the same time, only moderately 

dependent on the improvement prescription with &us2 l^us 2* = 0.82. On the basis of Eq. (48) this moment might 
well serve for a reliable m s determination, with a sufficiently careful interpretation of the theoretical uncertainty. 

The corresponding spin zero part is, per se, proportional to m 2 s (not counting non-perturbative, so-called "conden- 
sate" contributions) and thus could be considered as ideal for a measurement of m s . However, the behaviour of the 
perturbative series 

& = + 9.33a+110a2 + 1323 o2 



a^(12200 + 4 ll4 ±17.5£; 2 2]3 ) 



mi 



r „ (1. + 0.992 + 1.24 + 1.59 + 2.16) 



3 m 2 
2M. 



oTF|(7-0±2) (49) 



and 



~C = 1^.19 + 11.2a +126. a 2 + 289. al + 6.63 : 

4(2.7l4 1]4 + 7.76fc 2 214 )) 



3 m 2 

-—|(3.19 + 1.19 + 1.42 + 1.94 + 2.6 ± 1.31) 

2 M± 

3 m 2 

2^(10.3 ±2.6) (50) 
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shows a rapid growth of the coefficients. The series is not expected to provide an accurate prediction for the mass 
effects. 



V. SUMMARY 



Implications of the newly calculated a^n 2 terms together with the a^n^ terms for an improved extraction of a s 
from the r decay are presented. Arguments are presented in support of predictions for the remaining terms of order 
a\nf and a A s n® which are based on FAC or PMS optimization. The complete calculation will lead to a reduction of 
the theory uncertainty within the frameworks of FOPT or CIPT down to a negligible amount. However, an irreducible 
difference between the results from these two schemes of 5a s (M T ) « 0.02 corresponding to Sa s (Mz) ~ 0.002 persists 
even after inclusion of 0{a i s ) (and even 0(a 5 8 ) terms). Similar investigations, based on data up to higher energies (e. 
g. for fictitious heavy lepton of 3 GeV of for sum rules based on e + e _ data) would lead to significantly smaller errors. 

New contributions of orders 0(m 2 a A s n 2 ) and 0(a 4 n 2 ) to (axial)vector correlators relevant for the QCD description 

of the semileptonic r decay into hadrons are obtained. The moments i?( 0, °) and i?( 0,1 ) are evaluated separately for spin 
zero and spin one final states The results are tested against predictions of FAC/PMS optimization methods. Good 
agreement is found. This has motivated us to take the full of FAC/PMS predictions as the basis for a new extraction 
of a s and m s from r-decays with 0(a A s ) accuracy. Using 8 C p V — 0.207 we find 0.330 and 0.1192 for af OPT (M T ) and 
a^ OFT (Mz) respectively. In the framework of contour improved evaluation these values increase to 0.354 and 0.1219 
respectively 

In contrast to the massless result, the PT series contributing to the m 2 dependent part seem to be barely convergent 
and the additional higher order terms seemingly do not lead to any significant improvement of the theoretical accuracy 
in the determination of the strange quark mass from the r decays. 

A slightly more favorable pattern of convergence is observed for the moments of the spin one contribution separately. 
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[40] This is in agreement to the analysis of Ref. where it has been concluded that "the resummed values of a a from r 
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